Origin of Diffuse Scattering in Relaxor Ferroelectrics 
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^S) , High-pressure and variable temperature single crystal synchrotron X-ray measurements combined 

with first-principles based molecular dynamics simulations study diffuse scattering in the relaxor 
ferroelectric system PSN (PbSci/2Nbi/203). Constant temperature experiments show pressure 
induced transition to the relaxor phase at different temperatures characterized by butterfly and rod 
shaped diffuse scattering around the {hOO} and {hhO} Bragg spots, respectively. The simulations 
' reproduce the observed diffuse scattering features as well as their pressure-temperature behavior, 

and show that they arise from polarization correlations between chemically-ordered regions, which 
in previous simulations were shown to behave as polar nanoregions. Simulations also exhibit radial 
diffuse scattering (elongated towards and away from Q=(000)), that persists even in the paraelectric 
' phase, consistent with previous neutron experiments on (PbMgi/3Nb2/303) (PMN). 

C/3 ! 

Single crystal relaxors have huge electromechanical coupling, and show much promise for ultrasonic transducer 
applications a. They have broad frequency and temperature dependent dielectric maxima, which in the special case 
of relaxor- ferroelectrics 0,0 drops to a much lower value below the ferroelectric transition temperature (Tpg). The 
origin of the relaxor phase has been a topic of intense research for over a decade. From refractive index measurements, 
Burns et al. Q suggested that formation of polar clusters below a characteristic temperature, now called the Burns 
I temperature (T^) [J, gives rise to dielectric dispersion, which was confirmed experimentally Q. These clusters thought 

■ to be a few nanometers in size are called polar nanoregions (PNRs). 
C \ Recent X-ray and neutron experiments show characteristic shapes of diffuse scattering in the relaxor phase of 
^ ■ several lead based relaxors [1, 0, ll, 0] absent in their paraelectric or ferroelectric phases. The main observed feature 
I— — 1' is anisotropic diffuse scattering around the Bragg peaks along <110> directions. The ferroelectric phase shows weak 
streaks similar to BaTiOs and KNbOs, but rotated by 45-degrees. The paraelectric phase only shows radial diffuse 
scattering. 



X 



> ; Several hypotheses previously advanced to explain the characteristic shapes [6l Ufl LUl , invoke some type of 
\ artificial anisotropic features such as anisotropic strain, correlated atomic shifts, or domain boundaries to generate the 
■ experimentally observed anisotropic diffuse scattering features. Fitting the shapes of diffuse scattering features is not 
sufficient to uniquely determine the nature of the microstructural feature that give rise to them. Incorporating realistic 
00 polarization fluctuations via first-principles derived models provide a basis for clarifying the nature of the PNR in 
relaxors and experimental diffuse scattering observations provide a critical test for any theoretical model. In this 
Letter, we clarify the microstructural origin of diffuse scattering features in the relaxor- and the paraelectric- phases 
r\ , of PSN, and related materials, by combining single crystal X-ray diffraction experiments and first-principles-based 
molecular-dynamics simulations [l3| . Our result suggests new pathways to engineer novel materials with superior 
electro-mechanical properties. 

X-ray diffraction was measured on single crystal PSN at beamlines 11-ID-C (0.1077 A) and 16-BMD (0.436926 A) 
\ of the Advanced Photon Sources (Argonne National Laboratory). Single crystals of disordered PSN were grown from 
high-temperature solutions using a mixture of PbO and B2O3 as solvent [1^. A crystal with dimension of 70x70x20 
/xm^ and oriented along (001) was loaded into a diamond anvil cell (DAC) with Ne as the pressure medium. A ruby 
chip and a small grain of Au were also loaded for pressure determination (accuracy of 0.2 GPa). A MAR3450 image 
plate was used to record oscillation photographs. For low temperature measurements, the DAC was loaded in a He 
flown cryostat, in which temperature was measured with a thermocouple with an accuracy of ± 2 K. Additional 
experimental details about beamline 11-ID-C and 16-BMD can be found in Refs. [l^, To facilitate measurement 
of the diffuse scattering signal, the sample was rocked ± 6° in omega. Because of its (001) orientation in a cryostat, 
only the (hkO) indexed Bragg peaks could be observed. 

At ambient pressure, PSN enters the relaxor phase below the Burns temperature ( T^ ~ 650 K) and undergoes a 
ferroelectric transition at TpE ~ 365 K. In our experiments at 400 K (Fig. la), we observe the characteristic butterfly 
and rod shaped diffuse scattering. At 250 K (Fig. lb), below TFE^ the characteristic diffuse scattering diminishes 
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FIG. 1: Temperature dependence of diffuse scattering in (hkO) plane from X-ray diffraction at ambient pressure (red background) 
and simulations at P=18GPa (white background) (a) Characteristic butterfly and rod shaped diffuse scattering (enhanced 
intensity along [110] and [110] directions avoiding Q =(000)) is observed in the relaxor phase around {hOO} and {hhO} Bragg 
spots, respectively, in experiments (400 K) and simulations (200 K) (b) In the ferroelectric phase, the characteristic shapes 
diminish in experiments (250 K), similar to simulations (10 K) (c) In the paraelectric phase, simulations (320 K) only show 
radial diffuse scattering (d) Ice only shows the butterfly and rod shapes, while Idd shows only the radial patterns in the simulated 
relaxor phase. 



considerably. Diffuse streaks connecting the Bragg peaks grow weaker, yet persist in the ferroelectric phase, similar 
to those observed in BaTiOa and KNbOs (Ref. [17| and references therein). 

We observe pressure induced phase transitions at 300 K, where around 1.4GPa a reappearance of the relaxor phase 
characterized by diffuse scattering is seen; persisting up to l.SGPa, but disappearing at l.SGPa (Fig. 2b). Disappear- 
ance of diffuse scattering at higher-pressures indicate a ferroelectric to paraelectric phase transition consistent with 
dielectric- measurements . 

At 50 K (Fig. 2b) we also observe superlattice peaks at low pressures of P=0.6 GPa at (h-f 1/2 k+1/2 0) suggesting 
a lowering of symmetry, which reduce in intensity as pressure drives the system to the relaxor phase. Our experiments 
are on disordered PSN, hence the peaks are not due to ordering but due to possible octahedral rotations. In the 
relaxor phase, at 2.7 GPa, additional satellite peaks with diffuse wings emerge, which persist up to 10 GPa. Ordering 
PNRs with a range of closely spaced incommensurate modulation vectors will give similar satellite spots with diffuse 
wings. 

We performed molecular-dynamics simulations using a model Hamiltonian obtained by expanding the potential 
energy of PSN about a high symmetry perovskite reference structure and projecting onto the subspace of soft normal 
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modes, dominated by Pb displacements, including the ferroelectric instabilities 1^ 2^- Model parameters are fitted 
to first-principles density functional calculations [2l|, [23] ■ 

The model has chemically ordered regions (COR) with rock-salt ordering of "Sc+'^" and "Nb+^" 'B' site ions, 
embedded in a chemically disordered region (CDR) that has a random 'B' site configuration. This gives rise to 
a quenched "random" component to the local electric fields at the Pb-sites. Relative to equal volumes of CDR, 
polarization is enhanced in CORs, thereby acting as PNRs. The CORs are spatially fixed, because chemical order 
(disorder) is quenched, but their average polar orientations vary dynamically [2l|. The model allows homogeneous 
strain to fluctuate; inhomogeneous strain, higher-frequency contributions to lattice polarization, and oxygen octahedral 
tilting are ignored. 

In X-ray and neutron scattering experiments, the measured intensity is a Fourier transform of the ensemble averaged 
two point density-density correlation. Assuming point particles, the intensity is: /(Q) =< ^ e'*^ '^'''"'"^^ >=< 
I Si giQ-'^i^s > where i runs over the total number of Pb-atoms in the supercell (Pb scattering is expected to dominate 
our experiments as it has a disproportionately large X-ray scattering factor) and Q values commensurate with the 
supercell were chosen. Because our distribution of COR in our supercell is not entirely isotropic, we apply full cubic 
symmetry to the scattering intensity before comparing with experiments. 

Fig. 1 shows the temperature dependence of our computed diffuse pattern in the (hkO) plane for the ferroelectric, 
relaxor and the paraelectric phases at P=18 GPa. In the relaxor phase (Fig. la) we observe the characteristic 
butterfly and rod shaped diffuse scattering around (hOO) and (hhO) Bragg spots, which diminish in their intensity in 
the ferroelectric phase (Fig. lb) in excellent comparison to experiments. 

Pressure dependent changes in diffuse scattering (Fig. 2a) are similar to those induced by varying temperature. 
Again, our model captures the experimental pressure induced phase-transition at BOOK (Fi g. 2 b). Experimentally 
observed streaks, like in KNbOa have been shown to be due to hopping between equivalent sites [23|, and superstructure 
peaks at 50 K are possibly due to octahedral rotations, both of which are excluded from our model. But these atomic 
features do not cause the relaxor phase, which is characterized by the butterfly and rod shaped scattering, which our 
simple model is able to capture. 

As the system is driven towards the paraelectric phase, by increasing temperature or pressure, the characteristic 
butterfly and rod shapes vanish; however, a weak diffuse pattern that extend radially towards and away from the 
origin (Q=(000)) persists around all the Bragg peaks. Its iiitensity is larger in the direction away from the center 
than towards it, consistent with recent neutron experiments [9] but absent in X-ray due to their low intensity. 

We find no significant evolution of the diffraction patterns during our simulations, which access frequencies greater 
than lOGHz. This suggests only static correlations cause diffuse scattering, consistent with recent experiments on 
PMN [9|]. 

We find the contributions to the scattering intensity from correlations due to atoms in the different regions (COR 
and CDR). These are shown in Fig Id, where Ice and Idd are the intensities obtained by setting CDR and COR Pb 
displacements to zero respectively. While Ice shows characteristic rod and butterfly features, the radial diffuse pattern 
is absent. Idd only shows weak radial diffuse pattern. Thus, it is the correlation between COR atoms that give rise 
to the characteristic diffuse pattern. The CDR Pb-displacements are dominated by strong electric flelds flxed in a 
"quenched" distribution of "Sc" and "Nb" atoms in the chemically disordered matrix, and therefore radial scattering 
persists even in the paraelectric phase. This is consistent with experimental observations on PMN Note that 
our interpretation that radial diffuse scattering is associated with chemical disorder differs from the interpretation of 
Gehring et al. that it indicates chemical short range order 

To understand the origin of the shapes of the diffuse scattering we write the total intensity as a sum of Bragg and 
diffuse scattering, and expand the latter in powers of — ri — Roi, which being proportional to the local polarization 
is a small quantity. The Bragg term is j giQ (Roi-R-oj)^ aj^^]^ ^}jg lowest order diffuse scattering term is: 

Id^ff{G,<i) =< ^ |G + q|2(^, • a)(0 • a) cos(Q • (Rq, - Rq,)) > (1) 

where q=Q-G and G is the Bragg spot, a is the unit vector along G + q. The summation of pairs of atoms 
is equivalently re- written as one over inter-atomic distances R — Roi — Rqj , so that Eq. [T] becomes the Fourier 
transform of projections of real space displacement-displacement auto-correlation tensor C(R) along a: 

lMff{G,q) =<^(a-C(R)-a)|Q|2cos(Q-R) >, (2) 

R 
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FIG. 2: Pressure dependence of diffuse scattering in simulations at 180 K (white background) and experiments (red background), 
(a) Simulated diffuse scattering in the ferroelectric (P=14 GPa), relaxor (P=18 GPa) and the paraelectric (P=25 GPa) phases 
respectively, show characteristic relaxor-phase butterfly and rod shaped diffuse scatting and radial diffuse scattering in the 
pressure induced paraelectric phase (b) Experiments at 300 K show weak diffuse scattering in the ferroelectric phase (P=0 GPa) 
but strong butterfly and rod shaped scattering in the relaxor phase (P=1.4 GPa). In the paraelectric phase (P=1.8 GPa) the 
butterfly and rod shapes disappear, (c) At 50 K, we observe superlattice peaks in the ferroelectric phase (P=0.6 GPa). In the 
relaxor phase (P=2.7 GPa) on top of the butterfly and rod shapes we observe satellite spots. Further increase in pressure to 
10 GPa destroys the diffuse scattering, leading to a pressure induced relaxor to paraelectric phase transition. (The bright spot 
at the center of the image is due to the beam stopper.) 



where C(R), given, in symmetric form is 

Cq/3(R) = - y^(CiaCz±R;3 + ^i^CiiRa), (3) 

i 

and the "±" notation impHes summation over both signs. 

For small |q|, we approximate G + q by G and a by the unit vector do- Fig 3. shows the real space correlations 
in the relaxor phase (180K, 18GPa) with projections along Q!g=[100] and [110] summed over the 'z' direction as 
well as their Fourier-transforms. The diffuse scattering butterfly and rod shapes come from aG = [100] and q;g = [110] 
projections, respectively, of the COR-COR real-space correlations. The CDR-CDR regions have strong correlation 
along the direction perpendicular to the radial direction, leading to the radial diffuse scattering in the Fourier space. 
In 3D (not shown) the [110] (as well as [111]) real space projections appear as ellipsoids, while the [100] correlations 
appear as discs. 
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FIG. 3: Real-space correlations (Eq. [Sjl projected along [100] and [110] directions and their direct Fourier transforms, (a) and 
(d) show correlations between all 'Pb' atoms, which reproduce the experimentally observed butterfly- and rod- shaped diffuse 
scattering, (b) and (e) show the correlations between COR-COR regions which are extended in the real space along directions 
normal to the butterfly and rod directions. Their Fourier transform gives the characteristic diffuse scattering, (c) and (f) show 
the correlations between the CDR-CDR regions which when Fourier transformed give the radial patterns 



We conclude that PSN shows ferroelectric —+ relaxor ^ paraelectric phase transition with increase in temperature 
and/or pressure in our experiments. The relaxor phase in the experiments is characterized by butterfly and rod shaped 
diffuse scattering. 

Simulations further reveal that the anisotropic correlations from COR-COR 'Pb' atoms which only have orienta- 
tional degrees of freedom, give rise to the characteristic relaxor-phase diffuse scattering, rather than the effect of strain 
or artificial atomic shifts or domain boundaries. Polarization would be accompanied by local strain inhomogeneities, 
that would cause additional contribution to the diffuse scattering [ijj, but are a secondary effect. The radial diffuse 
scattering is identified as coming from local concentration fluctuations at the B-site, which reflect in the Pb displace- 
ments from corresponding fluctuations in local electric-fields on Pb. Coulomb energy is minimized when the local 
polarization aligns parallel to the local field. Our result suggests that by engineering the shapes and relative positions 
of the COR regions, one could possibly control the anisotropy in the COR-COR correlations, leading to the design of 
new relaxor materials with superior electromechanical properties. 

An increase in pressure reduces the ferroelectric (free energy) well depth and therefore the ferroelectric Pb- 
displacements in COR. This reduces PNR-PNR correlations (and not their sizes) and hence the butterfly and rod 
shaped diffuse scattering. In the ferroelectric phase, the whole system, including the disordered matrix, has a uniform 
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polarization along <111> directions. This enhancement of the local polarization reduces the effect of correlated Pb- 
displacements in COR-COR pairs, greatly diminishing the characteristic butterfly- and rod- shaped diffuse scattering 
features. With increase in temperature, thermal fluctuations of the polarization decrease the magnitude of the COR 
polarization, thereby decreasing the COR-COR correlations. Above T^, only the CDR-CDR radial contributions 
remain. 
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